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q \ We present a higher-twist production mechanism for r\ c and xj charmonium at large 

momentum fraction xp in pion-nucleon collisions. The higher-twist contribution is 
essentially independent of xp and therefore becomes dominant at large xp, where 
the leading-twist contribution falls off as (1 — xf) 3 . We show that the higher-twist 
mechanism produces longitudinally polarized Xi an d Xi- For the X2, this is clearly 
different from the leading- twist prediction of transverse polarization. For the xi> the 
^ . polarization of the leading- and higher-twist contributions is qualitatively similar. 
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1 Introduction 



Heavy quark-antiquark bound states (quarkonium) were among the first systems to 
be successfully studied using perturbative QCD. Much of the quarkonium decay data 
was adequately explained by the colour-singlet model [p], §, |J. In this model, the 
quarkonium states are taken to be non-relativistic colour-singlet QQ states, and decay 
amplitudes are written as convolutions of the quarkonium wavefunction with the 
lowest-order perturbative amplitudes for QQ annihilation. It could then be expected 
that quarkonium states are dominantly produced by the leading-twist mechanisms 
where the subprocess {e.g. gg — > J/ip + g) is related to a decay mechanism {e.g. 
J/ip — > ggg) by crossing @, || . However, the resulting predictions are contradicted by 
several experiments, from fixed-target [|, [7], || to the highest collider energies |], |10| . 

In this article, we shall concentrate on the production of charmonium at low trans- 
verse momentum P± and large longitudinal momentum fraction xp in fixed-target 
pion-nucleon collisions. We have earlier pointed out that the anomalous dependence 
of the production cross section on the mass number A of a nuclear target [|11[] suggests 
the importance of higher- twist, non-factorizable production mechanisms of the J/ip 



in this region f!2 |. In a recent work |]13 , we calculated the polarization of J/^'s 
produced in pion-nucleon reactions in the colour-singlet model at leading twist, in- 
cluding contributions from direct parton fusion into J/ip and from radiative decays 
of the Xi,2 charmonium states. We found a discrepancy between these leading-twist 
predictions and experimental data || [?[]. The purpose of this paper is to present 
a specific higher- twist mechanism for the production of C = +1 charmonia (the r\ c 
and xj) a ^ large Xp. We shall show that higher-twist mechanisms are likely to be 
dominant at large xp and that their characteristic features may explain some of the 
charmonium polarization data. 

The polarization of J / ip charmonium produced in pion-nucleon collisions has been 
measured by the Chicago-Iowa-Princeton collaboration M for 0.3 < xp < 1.0. They 
observe an abrupt change from unpolarized production to longitudinal polarization 
at xp ~ 0.95. This effect is reminiscent of the longitudinal polarization effect ||l4j, |15H 
observed in the production of dileptons in pion-nucleon collisions, which may be 
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explained by a higher-twist production mechanism first suggested by Berger and 
Brodsky 0, |17|, [L8|, [0|fg . 



The higher-twist charmonium production mechanism which we present here is 
closely related to the higher-twist Drell-Yan mechanism. We have chosen to study 
the C = +1 states because they couple to two gluons, which makes this the simplest 
case to calculate. We shall show that the xi,2 are produced longitudinally polarized 
at large xp. For the X2, this is strikingly different from the usual leading- twist gluon- 



fusion mechanism, gg — > X2-, which only produces transversely polarized X2 |2T 
For the Xi, on the other hand, the leading- and higher- twist contributions have similar 
polarization properties, so that the signature of the higher-twist mechanisms is not 
as clear as for the X2- 

Thus a measurement of the polarization of the X2 will provide important informa- 
tion about the production dynamics of charmonium. An observation of the predicted 
longitudinal polarization of the xi,2 at large xf would support the helicity conserva- 
tion rule between the pion and the large xp system |IB |, which is already supported 



by data for the J/ip and the Drell-Yan virtual photons. 

In Section 2 of this article we outline the calculation of the charmonium production 
cross sections. We illustrate their xf dependence in 7] c production, in which case 
compact analytical expressions can be derived. In Section 3, we present numerical 
results for the polarization of the Xi,2- in Section 4, we summarize our results and 
discuss the relation of this model to other models of charmonium production. 

2 Outline of the calculation 

The production of systems that carry a large momentum fraction xp requires a pertur- 
bative momentum exchange involving all the valence partons of the projectile |16|, |22 



Thus a correlation scale is present in the reaction. The cross sections become sup- 
pressed by powers of the ratio of this correlation scale and the hard momentum scale, 
e.g. f^/m 2 c in charm production from pions. These mechanisms are therefore of higher 
twist and are usually neglected. At large xf, however, inverse powers of (1 — xf) may 
compensate the suppression, and higher-twist mechanisms may become dominant 



3 



Here, we present a higher- twist mechanism for C = +1 charmonium production. 
Our mechanism is described by the Feynman diagram in Fig. 1. Note that whereas the 
mechanism of Ref. |TjJ produces both the higher-twist component of the cross section 
and the perturbative tail of the leading-twist component, our mechanism is purely 
higher twist. This is evident from the Feynman diagrams. As in the Drell-Yan case, 
we treat the interaction with the nuclear target simply as the annihilation of one of 
the pion's valence quarks with an on-shell quark from the target. This is the simplest 
higher-twist mechanism for charmonium production, but not necessarily the most 
important one. There may be a significant contribution from higher-order processes 
where a gluon from the target interacts with the pion, because at high energies and 
large Xp the parton distribution of the target is sampled at small ^target- Moreover, 
the simplifying assumption that the target just contributes one parton implies that 
the cross section is factorizable in terms of the target parton distribution. An analysis 
|T^ 1 of J/ip production data JTIJ suggests that there are important non-factorizable 
production mechanisms. 



We treat the pion bound-state effects as in Ref. |24|] and the charmonium binding 
as in the colour- singlet model |25fl . In other words, the amplitude for the reaction 
7r~ +u — ► d + 2S+l Lj is written as a convolution of three factors: the pion distribution 
amplitude 4>(z), the perturbative-QCD amplitude for ud + u — > d + cc, and the non- 
relativistic wavefunction $(P, q) of the charmonium state. Here, z is the light-cone 
momentum fraction carried by the u quark in the pion, P is the total four-momentum 
of the charmonium state, and 2q is the relative four-momentum of the charm quarks. 
The amplitude is 

A(n~ + u - d + 2S+1 Lj) = Ana ^ u c { Pd )r 2 l,i>^u c { Pu ) 



J zs(l — z)u 

where s = (p n +p u ) 2 , u = (p n —Pd) 2 , C is a colour index, and A MlAi2 is the truncated 
amplitude for gg — > 2S+1 Lj, which is the convolution of the quarkonium wavefunction 
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$(P, q) with the hard amplitude M1 u 2 for gg — > cc 
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as explained in Ref. PS| (an example is given in eq. (|TID below). The convolutions 
of the pion distribution amplitude with 0^ lfJj2 and its derivative with respect to q, 
evaluated at the non-relativistic limit q = 0, are 
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The singularities at 2 = 0, 1 from the gluon propagators are expected to be cancelled 
by the distribution amplitude <p(z), whereas the singularity at z = zq from the heavy 
quark propagator is regularized by the +ie prescription, which gives the amplitude a 
non-trivial phase [p6|1 . Below, we shall present results for the symmetric (asymptotic) 
24]] and two-humped [^7j distribution amplitudes^, 



6z(l-z), 
z(l -z) [26 



100z(l -z)\. 



(9) 
(10) 



1 The simplest model for the pion's valence state, <p(z) — S(z — 1/2), leads to rather non-realistic 
shapes of the cross sections since the +ie prescription fails to regularize the singularity of the quark 
propagator if the distribution amplitude does not vary smoothly over the pole. The same is actually 
true of the higher- twist contribution |16 to the Drell-Yan cross section, although the singularity is 
cancelled when one considers angular distributions. 
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We shall first study r/ c production with a symmetric pion distribution amplitude, 
in which case the analytic expressions are quite compact. For the i] c , the relation 



between v4 M1M2 and C M1M2 is [0 



A-ti 



R s (o) 



where M = 2m c is the mass of the charmonium state, and -Rs(O) is the value of the 
S-wave wavefunction at the spatial origin. The square of the nu — > dr) c amplitude, 
averaged over the spin and colour of the u quark, is 

32(4t™ s ) 4 |P 5 | 2 / 2 1 
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where L = ln[(s — M 2 )/(M 2 — u)]. The differential cross section is 
da 1 



l + Zo L) 2 +7T 2 ^]}, (12) 
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where s to t = (Pn + £>nucicon) 2 - In this and the following formulae, we use the light-cone 
momentum fraction x = (P° + P 3 ) / (p^ + p 3 ) instead of the longitudinal momentum 
fractional. There is a simple relation between these two: xp = x— (M 2 +Pf)/(xs tot ). 
The two contributions are from tt + u — > d + 2S+1 Lj and tt + d — > u + 2,s+1 Lj, 
respectively. 

In the limit M 2 ^> P 2 ^ with x finite, the cross section simplifies to 
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Now, taking the limit x — > 1 in ([Hp means taking the double limit 



M 2 > P 



1 with M 2 (l - x) > P 



2 



(15) 



i.e. the same limit as for the Drell-Yan process in Ref. fll6 |. Although the limit 



(15) is not quite reached in present-day pion-nucleon reactions [20], it illustrates the 
behaviour of the Pj_-integrated cross section at large x since it takes into account the 
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fact that as x is increased, the dominant region of the integration shifts to smaller 
P±. The cross section fll4|) simplifies to 



da 512Tr 2 aj\R s \ 2 rt 



fu/N WJ + h,N w) 



1 

p 



(16) 



dxdPl 27M 5 s to t 
The integral / dP]_ over the perturbative region where P±/{1 — x) ^> Aq CD brings 
a logarithmic factor. Apart from this logarithm, the higher-twist contribution is 
suppressed by f 2 /M 2 with respect to the leading- twist, leading-order contribution 
from the subprocess gg — > r\ c |J, which is 

dP ±lZim = -Tils— - fg/MfgM—-). (17) 



o 



'dxdP 2 3M 3 s tot x 9 9 xs tot 
From ( |16|) we also see that da/dx remains constant as x — > 1, up to the logarithmic 
factor from the P± integration. The leading-twist contribution, on the other hand, 
falls off according to the power behaviour of the pion's gluon distribution, i.e. (1 — a;) 3 
from the spectator counting rulesQ ||28|| . At very large x, this will compensate the 
suppression by f 2 /M 2 , and the higher- twist component will become dominant. 

3 The polarization of charmonium 

We now consider the production of the Xi,2 states within our higher-twist model. The 
cross sections for absolute values of the helicity |A| = 0, 1, 2 can be extracted by using 
the covariant polarization sums given in e.g. Ref. |[32|| . However, the analytical ex- 



pressions that are obtained by using a symbolic manipulation program such as Reduce 



33] are too lengthy to be given here (some of them contain over a thousand terms) 
or even to be used in numerical integration routines. Instead, we have evaluated the 
double-differential polarized cross sections da/dxdP]_ (|A| = 0,1,2) numerically at a 
characteristic value of the transverse momentum, P]_ = M 2 (l — x). 
In Fig. 2 we plot the ratio 

Ea=±i da/dxdPl (A) , . 



R(l) 



da/dxdPl (A = 0) 



Experimental determinations of the large x behaviour of the gluon distribution [^9], ||o|, |3l| arc 
based on fitting a purely leading- twist model to Drell-Yan, prompt-photon and/or charmonium data. 
The resulting powers of (1 — x) therefore partly reflect a hardening of the cross section due to the 
higher- twist component. 
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of the polarized cross sections for the Xi as a function of xp, using the symmetric 
and two-humped distribution amplitudes. The ratio of the polarized leading-twist 



cross sections da/dx [13| is also shown. In Fig. 3 we plot the ratios R(l) and R(2) 
for the \i- The leading- twist, leading-order gluon-fusion mechanism gg —> xi gives 
da/dx (0, ±1) = 0, up to corrections of the order of 15 % from transverse momentum 
smearing of the pion fli"3 |. 

As xp — > 1, the contribution from our mechanism to both the Xi an d X2 cross 
sections becomes longitudinally polarized (A = 0). For the xi-, this is in striking 
contrast to the leading- twist prediction. For the xii however, the polarization of the 
higher-twist component is qualitatively similar to that of the leading-twist component, 
so that polarization cannot be used to discriminate between leading- and higher-twist 
mechanisms as for the X2- 



4 Discussion 

When a charmonium state is produced at a very large momentum fraction xp, all of 
the projectile hadron's momentum must be transferred to the charm quarks. This 
requires a higher-twist mechanism, where all the valence quarks of the projectile 
couple perturbatively to the charm quarks. The higher-twist nature of these processes 
is due to the soft correlation scale, e.g. f n , from the integration over the transverse 
momentum distribution of the valence state. In this article, we have presented a 
higher-twist mechanism for the production of the t] c and x J charmonia in pion-nucleon 
collisions. The contribution from this mechanism dominates over the leading-twist 
contribution at large xp. As seen from eq. flT5| ) for the r] c , the higher-twist cross 
section is almost flat as x — * 1, whereas the leading-twist cross section falls off as 
(1 -x) 3 . 

In our model, both the Xi an d Xi are produced longitudinally polarized at large 
xp. For the X2, this result is in striking contrast to the leading- twist prediction of 
purely transverse polarization; for the Xii the polarization of the leading- and higher- 
twist contributions is qualitatively similar. A direct measurement of the polarization 
of the Xi would give important information about the production mechanisms of 
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charmoniumQ. 

The mechanism considered in this paper is the simplest but not necessarily the 
most important of many possible higher-twist production mechanisms. In mechanisms 
such as this and the Drell-Yan mechanism of Ref. JTBI , where the valence state of 
the pion interacts with a single parton from the target to produce a large xp system 
in the final state, the parton distribution of the target is at high energies sampled 
at small values of x target . Due to the dominance of gluon distributions over quark 
distributions at small ^target, there may be significant contributions from higher-order 
(in ctg) mechanisms where the target contributes a gluon instead of a quark. 

Mechanisms where the target contributes a single parton are factorizable in terms 
of the parton distribution of the target. Thus they would not bring about the vi- 



olation of factorization seen in the reaction n + A — > J ftp + X \\12\\ . Instead, the 
nuclear number dependence of the cross section would be consistent with shadowing, 
as is experimentally the case in the Drell-Yan reaction. More general higher-twist 
mechanisms outlined in Ref. [23], where a spectator valence quark interacts softly 
with the target, would give factorization- violating cross sections. 

In this article, we have compared the leading- and higher-twist components of 
the charmonium production cross sections within the colour-singlet model. There 
is also a potentially important contribution to the production of the xj states from 
mechanisms which involve non-perturbative transitions between colour-singlet and 
colour-octet QQ states p3 |. These mechanisms are necessary in a rigorous calcu- 
lation of the xj decay widths |56fl , which in the colour-singlet model are infrared 
divergent at 0(a^). In the rigorous analysis, there is a contribution proportional to 
a 1 < Xj\^s( 3 Si)\xj >, where the infrared singularity has been absorbed into the 
non-perturbative matrix element < Xj\^^Si)\xj >■ In Xj production, there is 
a corresponding contribution from the non-perturbative production matrix element 
< 0|Og J ( 3 S'i)|0 >. The magnitude of this matrix element can be measured in B- 



meson decays into xj charmonium [P? , 58fl . In our mechanism, where the xj states 

3 The decay angular distributions of P-wave charmonia have actually been measured recently 
pifl . However, no separation was made between the xi an d X2 i n the angular distribution analysis, 
and due to limited statistics, the errors are still quite large. With the present statistics, it is also 
impossible to determine the polarization as a function of the momentum fraction xf- 
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couple to two gluons, there is no infrared singularity that would signal the break- 
down of the colour-singlet model. Nevertheless, the octet contribution is of the same 
order in the perturbative and non-relativistic expansions as the singlet contribution. 
To estimate the importance of colour-octet mechanisms in xj production at higher 
twist, one should calculate the cross section for the process 7r~ + u — > 3 Si + d where 
the QQ( 3 Si) heavy quark-antiquark pair is in a colour-octet state. We postpone the 
analysis of these contributions to future work. 

In spite of these reservations, the present calculation illustrates the perturbative 
dynamics behind the hadron helicity conservation at large Xp, which has been ob- 
served in the Drell-Yan process |TJ], [T5| and in J/ip production p|. 

Finally, we note that the present mechanism and other similar mechanisms where 
the charmonium state couples to virtual gluons allow the production of xi & t the 
same order in a s as \2- in the leading twist case, the production of \i m the leading- 
order subprocess of on-shell gluon fusion, gg — > xi, is forbidden by the Landau- Yang 
theorem P9| , and consequently the ratio of the total Xi an d X2 production cross 



sections is predicted to be small. The experimental data shows equal cross sections M. 
If higher-twist production mechanisms are important down to low xp, as the failure 
of the leading-twist predictions of charmonium polarization |13| seems to indicate, 
they could explain the total cross section anomaly. 
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FIGURE CAPTIONS 



Figure 1. One of the two Feynman diagrams that describe the higher- twist 
charmonium production mechanism considered in this article. The four-momenta 
flow from left to right. The other diagram is obtained by crossing the gluon lines. 

Figure 2. The ratio R(l) = J2x=±idcr(X)/da(0) of the polarized xi production 
cross sections plotted as a function of xp. The higher-twist contributions have been 
evaluated at a characteristic value of the transverse momentum, P]_ = M 2 (l — x), 
using the symmetric and two-humped pion distribution amplitudes. The leading-twist 
contribution |Tj| has been integrated over P±. 

Figure 3. The ratios R(l) = J2 x=±1 da(X)/da(0) and R(2) = Ex=±2dcr(X)/da(0) 
of the polarized \2 production cross sections, plotted as a function of Xp. The higher- 
twist contributions, which are shown here, have been evaluated at a characteristic 
value of the transverse momentum, P\ = M 2 (l — x), using the symmetric and two- 
humped pion distribution amplitudes. The leading-twist prediction, with transverse 
momentum smearing and higher-order corrections neglected, is da(0, ±1) = 0. 
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